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Abstract. The interpretation of the Landau gauge lattice gluon propagator as 
a massive type bosonic propagator is investigated. Three different scenarios are 
discussed: i) an infrared constant gluon mass; ii) an ultraviolet constant gluon mass; 
iii) a momentum dependent mass. We find that the infrared data can be associated 
with a massive propagator up to momenta ~ 500 MeV, with a constant gluon mass of 
723(11) MeV, if one excludes the zero momentum gluon propagator from the analysis, 
or 648(7) MeV, if the zero momentum gluon propagator is included in the data sets. 
The ultraviolet lattice data is not compatible with a massive type propagator with a 
constant mass. The scenario of a momentum dependent gluon mass gives a decreasing 
mass with the momentum, which vanishes in the deep ultraviolet region. Furthermore, 
we show that the functional forms used to describe the decoupling like solution of 
the Dyson-Schwinger equations are compatible with the lattice data with similar mass 
scales. 
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1. Introduction and Motivation 

The lagrangian for pure SU(3) Yang-Mills theory does not include a mass scale. At the 
classical level, conformal invariance of pure gauge theories is the expression of this lack 
of scale. However, the corresponding quantum theory gets a mass, let us say Aqcd, 
from the loop contributions via dimensional transmutation. 

At the level of the lagrangian a gluon mass term is forbidden by gauge invariance 
and, as long as the color symmetry is unbroken, the gluon is supposed to be massless. 
Certainly, in what concerns the perturbative solution of QCD, within the framework 
of the Faddeev-Popov quantization procedure [1], the gluon is massless. However, as 
discussed in [2], a dynamical generated mass which is a function of the momentum is 
allowed if one goes beyond perturbation theory. 

From the theoretical point of view, a non-vanishing gluon mass is welcome to 
regularize infrared divergences and to solve some problems related with unitarity. On 
the other side, diffractive phenomena [3] and inclusive radiative decays of J/ip and T [1] 
suggest a massive gluon. Moreover, lattice simulations suggest an infrared gluon hard 
mass of ~ 600 MeV [5] and an ultraviolet mass Mg ~ 1.0 GeV [191 [20]. Phenomenology 
favors a gluon mass between ~ 0.500 GeV and ~ 1.2 GeV depending on how the mass 
is defined - see table 15 in [1]. Furthermore, a dynamically generated gluon mass can 
be related with the presence of the {A^) gluon condensate, which is associated with the 
non-perturbative sector of QCD and whose role has been investigated by several authors 
- see, for example, P, [3, [8] and references their in. 

The idea of a gluon mass was explored by different authors [9l [TOl [TTl [121 [ISl El 
[T5l [T6] . Starting from the Dyson-Schwinger equations (DSE) for the gluon and ghost 
propagators, after a suitable truncation scheme, the equations were solved and the 
transverse part of the gluon propagator described by a massive type propagator, i.e. 

with a a momentum dependent gluon mass M^(g^), called below running mass gluon, 
and a running dressing function Z{q'^). Typically, the numerical solution for the 
propagator is fitted to a functional form M'^{q^) copied from the solution discussed 
in [2]. The solutions of the DSE give a M{q^) which takes its largest value at zero 
momentum, where M(0) ~ 600 MeV, and vanishes for q ^ Aqcd- In this way, the 
usual perturbative propagator is recovered at high momentum. 

The Dyson-Schwinger studies referred to in the above paragraph rely on the 
Faddeev-Popov quantization method, whose validity for investigating non-perturbative 
effects in QCD is under debate - see, for example, pTl [T8] and references therein. 
However, a running gluon mass was also found within the studies of the non-perturbative 
quantization of Yang-Mills theories. In particular, in the framework of the so-called 
refined Gribov-Zwanziger action [18], a tree level propagator was computed suggesting 
a functional form for M(g^) which is close to the original proposal of Cornwall [2]. 

Besides the solutions of the Dyson-Schwinger equations, lattice simulations also 
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provide support for a non- vanishing gluon mass, see for example [SI [ISl EDI EI]. The 
precise value for M(g^) depends on how the gluon propagator is modeled. For example, 
in lattice QCD or Dyson-Schwinger calculations, a running mass is computed fitting 
the propagator to a given functional form for M(g^). In this work we use lattice QCD 
simulations to investigate M^(g^). Besides checking the compatibility of the lattice data 
with theoretical predictions for M^(g^), a first attempt is made to compute the running 
gluon mass directly from lattice simulations. 

The paper is organized as follows. In section [2] we describe the lattice setup, the cuts 
performed to produce a unique curve for the gluon propagator and the renormalization 
procedure. The gluon mass is investigated in section [3] considering three different 
scenarios. Before considering the running mass, we consider constant mass ansatze to 
fit the propagator for the infrared (section I3.1.ip and the ultraviolet momenta (section 
I3.1.2p . We find that the gluon propagator cannot be fully fitted with a constant gluon 
mass, and thus we proceed with fitting the gluon propagator with a running mass gluon 
and a running gluon dressing function (13. 2p . Finally, in section H] we resume the results 
of section [3] and comment on its interpretation. 



2. Definitions and Lattice Setup 

In this paper we consider four dimensional SU(3) Wilson action lattice simulations at 
P = 6.0. For this (3 value the lattice spacing, measured from the string tension [22], is 
given by a = 0.1016(25) fm or, equivalently, = 1.943 GeV. 

The gauge configurations were generated with the MILC code [23]. Each 
configuration, sampled with the Wilson action, was rotated to the Landau gauge, as 
described in [20], using for gauge fixing the overrelaxation algorithm. The gauge fixing 
process requires a maximization of a given functional over the gauge orbit of each 
configuration. The maximization process was stopped when the lattice tetra-divergence 
[20] . averaged per site, become smaller than 10"^'^. The set of lattices considered in the 
present work are summarized in table [H 

In the Landau gauge, the gluon propagator is given by 

D;W) = (A'^mt) = 6^' (5,. - D{q') ■ (2) 

latin letters stand for color indices and greek letters for space-time indices. The 
momentum space gluon field v4^(g) definition and how to compute the form factor D[q^) 
are described in [2^ and will not be repeated here. 

For the continuum momentum we take the standard definition 

= ^ sin (^-^ , n = 0,...,Lf,-l, (3) 

where a is the lattice spacing and the number of lattice points in direction fi. 

In the following only renormalized data will be considered. The renormalization 
was performed fitting, for each lattice simulation, the bare lattice propagator to the 
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Figure 1. Rcnormalizcd gluon propagator for all the lattices described in table [T] 
The propagator (vertical axis) is given in GeV~^ and the momenta (horizontal axis) 
in GeV. 



one-loop inspired result 



(4) 



where 7 = 13/22 is the gluon anomalous dimension. For the fits to equation (jl]), the 
largest momentum range, starting at qmin and going up to ~ 5 GeV, with the /d-of. 
closer to unity was used - see table [2l From the fits we extracted the constants K and 
A, which enabled the computation of the renormalization constant Zr via 

D{q^) = ZnDLaM^) , (5) 

after requiring the renormalized propagator to be given by 

1 



(6) 



As renormalization scale it was used = 3 GeV. 

All the simulations were performed on an hypercubic lattice which breaks rotational 
invariance. In order to reduce lattice spacing effects, for each lattice and for momenta 
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L 32 48 64 80 

L(fm) 3.2 4.9 6.5 8.1 

# Confs. 126 104 120 50 

Table 1. Lattice setup - the lattices considered are symmetric hypercubes. For 
conversion into physical units the lattice spacing we used a = 0.1016(25) fm, or 
= 1.943 GeV, as computed from the string tension |22j . 



L 


32 


48 


64 


80 


Qmin Qmax (GgV) 


2.81 - 5.08 


2.49 - 5.02 


1.51 - 5.14 


1.52 - 5.05 


xVd.of. 


0.91 


0.97 


0.89 


1.06 


Zr (GeV"2) 


0.149(21) 


0.149(18) 


0.1477(38) 


0.1478(53) 



Table 2. Fits of the conic cut propagator data to equation and renormalization 
constants as function of the the lattice volume, qmin and Qmax stand for the lowest 
and highest momenta, respectively, used in the fit to get Z^. 



q > 1 GeV the conic cut [TH] was applied. For momenta below 1 GeV, all the data 
points were considered. In this way, we hope to have a good description of the infrared 
region. The renormalization procedure, as described above, was performed separately 
for each lattice propagator. In all cases the fit to equation was smooth and the 
corresponding x^/d.o.f. ~ 1. The renormalized gluon propagator, after the cuts just 
described, is reported in figure [H As can be seen in figure [H finite volume effects are 
observed in our data. 

In the following, we will consider two different analysis of the lattice data: (i) a 
separate analysis for each lattice volume; (ii) the different lattices can be combined, 
after removing finite volume effects, to produce a propagator with a larger density of 
points in the momenta axis. If the different volumes allow for an evaluation of the finite 
volume effects, the combined data, with its larger density of points on the g-axis, will 
reduce the statistical error on the final result. 

Let us elaborate more on the combined data set. From figure [1] it is clear that 
the lattice propagator differ by more than one standard deviation as a function of 
the physical volume. This effect is clearly seen in the infrared region. To reduce the 
finite volume effects, we take the propagator from the largest volume considered here, 
corresponding to an 80"* lattice, as a reference and compare all the remain propagators 
to it. To reduce the finite volume effects, the infrared data of the smaller lattices was 
removed if the propagator was not compatible, within one standard deviation, with the 
80^ propagator. Due to the infrared cut, from the 64^ propagator only data with q > 425 
MeV was considered, from 48^ only q > 671 MeV data was considered and from 32^ 
only data with q > 848 MeV was included. 

The cuts just described do not remove all lattice spacing effects. Indeed, even 
after performing all the cuts one can observe data points where, for the same the 
propagator differ by more than one standard deviation. The difference is due to the 
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Renormalized Gluon Propagator - |l = 3 GeV 

P = 6.0 - combined 32^ 48"', 64"", 80* data 
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q [GcV] 



Figure 2. Renormalized gluon propagator D{q^). 

Renormalized Gluon Data - H = 3 GeV 

P = 6.0 - combined 32*. 48^ 64"*, 80* data 
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Figure 3. Renormalized gluon data for q^D{q^). 



violation of rotational invariance. Therefore, to minimize such type of effects, for each 
lattice volume and for the same coming from different g^, if the different estimates of 
the propagator don't agree within one standard deviation, one of the points is excluded. 
For example, for the 80^ lattice for momentum q = 457 MeV there are two estimates for 
the gluon propagator, D{q^) = 6.563(95) GeV~^ associated with the lattice direction 
n = (3,0,0,0) and D{q^) = 6.416(71) GeV"^ associated with the lattice direction 
77, = (2, 2, 1, 0), coming from different types of momenta. The first value is clearly above 
all the data points and it was not considered in the data sets. In this way, the surviving 
points will produce a unique curve for D{q^). 

The combined lattice data for gluon propagator, after performing all the cuts, is 
shown in figure [2l The corresponding q^D{q^) function is reported in figure O Figures 
E] and |3] suggest that both finite volume and finite lattice spacings have been removed 
from our data. 
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The largest lattice volume included in our simulation is (8.1 fm)^. Lattice gluon 
propagators computed with larger volumes, but with a lattice spacing about twice the 
spacing considered in the present work, were reported in [21]. The two sets of data 
were compared in [8]. We will not repeat this exercise but will resume the outcome of 
the comparison. The propagator used in this work and those of [21] are essentially the 
same. Indeed, they are compatible within one standard deviation for q > 200 MeV and 
show a small difference for smaller momenta. For q < 200 MeV the propagator of [23] 
is about 10% smaller than the data in figure [21 As seen in figure [H our simulations 
have a limited access to momenta below 200 MeV and, in this way, we expect that the 
impact of the finite volume effects on M^(g^) computed from the combined data set 
is well below 10% factor. Furthermore, the separate analysis of the each volume will 
provide an estimate of the finite volume effects on the running mass. 

Another source of systematics are Gribov copies, i.e. configurations which satisfy 
the Landau gauge condition but are related by finite gauge transformations. This 
is a difficult and computational very demanding problem for the lattice practitioner. 
However, the known SU(3) lattice simulations show that Gribov copies do not change 
significantly the gluon propagator, i.e. that the effect due to the copies are, typically, 
within the statistical error; see, for example, [20]. Therefore, in this work we will ignore 
possible effects due to the Gribov copies. 

3. The Gluon Mass 

Our goal is to obtain the gluon mass M. However, its computation, or the calculation 
of M^, is not independent of the estimation of the dressing function Z. In the next 
sections we look at the gluon mass as given by the gluon propagator, i.e. by equation 
([I]), and explore different definitions. First we consider constant mass ansatze to fit the 
propagator for the infrared and the ultraviolet momenta. However, the gluon propagator 
cannot be fitted over all momenta with a constant gluon mass. Then, we proceed fitting 
the gluon propagator with a momentum dependent mass, called below running mass 
gluon, and a running gluon dressing function. 

3.1. Constant Gluon Mass 

We start our analysis assuming a constant gluon mass, i.e. assuming that in equation 
M and Z are constants. Therefore, in this section we take the gluon propagator as 
being described by 

in a certain momentum range. 

3.1.1. Fitting of a Constant Infrared Cluon Mass To check if the infrared gluon 
propagator can be described by such type of model, the lattice data was fitted to 
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Figure 4. Results of fitting the infrared lattice propagators to equation ([7]) in the 
range [0,qmax] and for x^/d.o.f. < 2.0. In figure full symbols refers to fits including 
the zero momentum propagator, while open symbols mean that D{0) was removed 
from the data sets. 



equation ((Tj) in the momentum range [0,gmax]- The results are plotted in figure HI 
where full symbols mean that the zero momentum propagator was included in the fits, 
while open symbols mean that -D(O) was removed from the data sets. 

Figure H] shows a and Z that are, within one standard deviation, independent 
of the fitting range, i.e. of qmax- Furthermore, requiring a x^/d-O.f. < 2.0 means that 
the infrared propagator can be described by (JTj) for momenta up to g ~ 500 MeV, if 
one ignores the lattice estimate for -D(O), and up to g ~ 430 MeV, if the lattice -D(O) is 
included in the data set. More, only for volumes of ~ (6.5 fm)^ or larger, the infrared 
lattice propagator can be described by a massive type propagator if -D(O) is taken into 
account in the data sets. 

Values for Z, M and qmax for each data set are reported in table El The table 
includes only the results of the fits whose /d.of. is closer to unit. If one includes -D(O) 
in the data, Z and M seem to be essentially independent of the volume. However, if 
one removes the zero momentum propagator from the data set, then one can observe a 
slightly volume dependence with M and Z increasing with V . In both cases discussed, 
i.e. including of not -D(O), M and Z computed from the combined data is just below 
the corresponding values obtained from the two largest two volumes. 

For the fits where -D(O) is excluded, a linear extrapolation of M and Z with 1/L 
towards the infinite volume can be performed. Only after disregarding the 64^ results 
one is able to obtain excellent extrapolations. Note that the numbers reported for the 
fits to the 64^^ lattice propagator are off the main trend observed in the remaining fits. 
The linear extrapolation to the infinite volume give Z = 4.49(10) with /d.o.f. = 0.38 
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qmax (MeV) 


660 


508 


503 


505 


550 


z 


2.963(62) 


3.53(10) 


4.035(89) 


3.85(11) 


3.594(68) 


M (MeV) 


490(11) 


578(12) 


651(10) 


626(13) 


597(9) 


xVd.of. 


1.42 


0.99 


0.39 


0.95 


1.03 



Table 3. Results of fitting equation The values of M and Z reported are for the 
fits whose /d.o.f is closer to unit for each data set. 



and M = 723(11) MeV for a /d.o.f. = 0.70. To these extrapolated values corresponds 
a 

D{0)v=oo = 8.58(26) GeV~^ (8) 

in excellent agreement with the largest volume simulated here, -D(0)^^^g £^^^^4 = 
8.86(45) GeV^^, and with the infinite volume extrapolation computed in [8], -D(O) = 
8.3(5) GeV-2. Further, the value (E]) agrees with the -D(O) computed using large volumes 
and reported in [8]. 

Our conclusion being that the fits show that the infrared lattice propagator is 
well described by a massive type propagator from momenta up to ~ 500 MeV with 
an effective gluon mass around 650 — 700 MeV. From table [H it follows that if -D(O) 
is included in the data set, then a massive type propagator describes the lattice data 
up to moment ~ 500 MeV with Z = 4.044(78) and M = 648(7) MeV given by the 
average values of the 64^ and 80^ figures; errors computed assuming gaussian error 
propagation. On the other hand, if one excludes -D(O) from the data sets, again a 
massive type propagator describes well the infrared lattice data up to g ~ 500 MeV 
with Z = 4.49(10) and an effective gluon mass M = 723(11) MeV. 

3.1.2. Fitting a Constant Ultraviolet Gluon Mass The same reasoning applied to 
infrared can be used to investigate the high momenta region. However, for the high 
momenta the fits to ([7]) give a negative M^, with depending strongly on the fitting 
range. We take this result as an indication that the ultraviolet is not described by such 
a propagator. 

Our discussion of the high momentum region is not in contradiction with the results 
of [T9I [20] , where an ultraviolet gluon mass of ~ 1 GeV was claimed. In [191 EO] an 
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ultraviolet regulator was used and the lattice data surviving the conic cut fitted to 



iln(g2 + M2)(g-2 + M- 



(9) 



g2 + M2 

where M is the gluon mass. Notice that the positive gluon mass in the numerator of eq. 
([9]) is equivalent to a negative mass in the denominator of eq. ((71), and thus a negative 
ultraviolet mass, in the sense of eq. ([Tj), is not in contradistinction with perturbative 
QCD. 



3.2. Momentum Dependent Gluon Mass 

The lattice gluon propagator can be described by a massive type propagator in the 
infrared region. For ultraviolet momenta, the propagator follows closely the 1-loop 
perturbative prediction and, in this sense, one can claim that for the high momenta 
region the effective gluon mass vanishes. Therefore, a way of recovering this two results 
is through the introduction of an effective gluon mass which is a function of the gluon 
momenta. So let us assume that D^q^) is given by equation ([1]) and that M{q^) and 
the dressing function Z{q^) that are functions of the momentum. In the following 
we will refer to M{q^) and Z{q^) as the running mass and running dressing function, 
respectively. 

Our first try to compute M{q^) and Z{q^) is to reproduce the analysis of the two 
previous sections. After dividing the g-axis into small window momenta and, in each 
momenta window, fit the propagator assuming Z and M constant, the outcome is a 
decreasing M^{q'^) with and, for momenta above ~ 1 GeV, M'^{q^) becomes negative. 
Given that the computation of Z and are not independent, a negative mass squared 
is, certainly, due to a will conditioning of at least one of the functions. Further, the 
problem of a negative mass squared can be cured, for example, by changing the definition 
of Z. This requires some model building and, instead, we will proceed fitting the lattice 
data to well known functional formulas. 

Although, the outcome of this procedure to compute M'^{q^) and Z^q^) will not be 
reported here, we would like to call the reader attention to some results found when 
analyzed the running mass and dressing functions obtained as described in the previous 
paragraph. 

In what concerns M'^{q^), we find that in the infrared region it agrees well with the 
results reported in section [3.1.11 Further, M'^{q'^) follows a behavior in the infrared 
region and the lattice data can be fitted assuming M^i^q"^) ~ g^lng^ as suggest in [2]. 

The running gluon dressing function Z[(f'^ decreases linearly from g = up to 1 
GeV. In the ultraviolet region, ^'(g^) agrees well with the perturbative QCD prediction. 
We observed that the gluon dressing function is well described by the functional formula 

^(?') = 71 2TrT ' (10) 

^ ' [A + ln(g2 + mg)]^ ^ ' 
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Renormalized Inverse Gluon Propagator - |l = 3 GeV 

P = 6.0 - combined 32"", 48*. 64"', 80* dala 




q- [GeV-] 



Figure 5. Inverse gluon propagator for the combined data set. Note the "almost' 
linear behaviour with of 1/D{q^). 



where 7 = 13/22 is the anomalous gluon dimension and tUq ~ 1.6 GeV^. Note that in 
equation (fTO|) . plays a role similar to in In [191 120], the authors measured a 
M = 1.0(1) GeV. The fits give essentially the same value, i.e. an mo ~ 1.3 GeV. 

In figure |5] the inverse of the propagator is plotted, for the combined data set, 
against g^. The figure shows that 1/D[q^) follows essentially the perturbative behavior 
at large momenta, i.e. 1/D{q'^) oc g^, and at smaller momenta the data points deviate 
from a pure quadratic behavior. Not surprisingly, the observed ultraviolet confirms that 
D{q^) is well described by equation (j4]). 



3.2.1. M'^{q'^) from a Perturbative Inspired Z^q^): An effective gluon mass can be 
defined extending the perturbative behavior, observed in the high momentum region, 
towards the infrared region. Indeed, if one assumes that 



Z{q' 



log 1 



A2 



:iii 



is valid over all momentum range, which ensures that the gluon dressing function is 
finite everywhere, z and A can be measured fitting the ultraviolet lattice gluon data to 



D{q' 



log (1 + 1. 



q^ 



in a given range of momenta [qmin, Qmax] 
taking the lattice propagator as 
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91 

A2 



g2 + M2(g2 



(12) 

The running gluon mass is then defined by 



(13) 



for all q. With this definition, called below Mpertiq^)-, one ensures that the 
usual perturbative propagator is recovered for high momenta or, equivalently, that 
Mpertil'^) for sufficently high q. As described below, the gluon becomes massless 
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1.783(37) 


0.812(47) 


1.33 


68 




3.5 


1.904(63) 


0.661(72) 


1.13 


54 
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1.92(12) 


0.64(14) 


1.16 


39 



Table 4. Results of fitting equation for g G [qmin, Qmax] for (jmax = 5 GeV. 



for g ~ 1 GeV, which a surprisingly low mass scale. Of course, the exact momentum 
scale at which the gluon becomes massless is connected with the parameterization of 
the gluon dressing function Z{q^). 

The results of fitting the lattice data to equation f lT2|) . for q G [qmin,Qmax] and 
Qmax = 5 GeV, are summarized in table ID The reason to exclude momenta above 5 GeV 
is to avoid remaining lattice artifacts which where not removed by the renormalization 
procedure. In what concerns qmin, it should belong to a region where the perturbative 
behavior is recovered. The renormalization procedure described previously, see section 
|2] and table [2l suggests that g^m = 3 GeV belongs already to the perturbative region. 
The two remaining g^m values, i.e. 3.5 GeV and 4 GeV, where used to check for the 
stability of the fits. 

The results reported in table H] show that, for all the data sets, equation (fT2ll 
describes well the lattice propagators in all momenta windows considered. Further, 
within each data set, z and A computed for the different qmin are compatible within 
one standard deviations. There are a few exceptions, where the z and/or A become 
compatible with the figures of the fits having the larger number of d.o.f. within less than 
two standard deviations, namely the fits to the 32^ data with the smallest momentum 
range, the A from fits to the 48^ propagator with the smallest momentum range, the fits 
to 64"^ data with qmin = 3.5 GeV, the A from fits to the 80"^ propagator with qmin = 3.5 
GeV, the fits to the combined data set with qmm = 3.5 GeV. Moreover, comparing the 
fits with the larger number of degrees of freedom, it follows that z and A, for all data 
sets are, compatible within one standard deviation. Given the good stability of the fits, 
in the following, we will the results from the fits with the largest number of degrees of 
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Figure 6. Mpg^((g^) for the various data sets. Only statistical errors from the 
propagators are take into account. Statistical errors were computed assuming gaussian 
error propagation. 

freedom, i.e. for g^m = 3 GeV, to compute M^(g^). 

Figure |6] shows Mpert{<f) for all data sets. The errors on Mpg^^(g^) include only the 
contribution form the lattice propagator and were computed assuming gaussian error 
propagation. Mpertiq"^) starts at value ~ 1.2 GeV for q = 153 MeV and decreases when 
increases. For all the data sets and all fits reported in table HI Mp^rtiq'^) vanishes for 
momenta around 1 GeV. For the last point included in figure Ej the error bar is of the 
same order of magnitude as the central value. For momenta above 1 GeV, not shown 
in figure [6l M^^^.^^q'^) starts to fluctuate and the error bars increase. 

The data for the different lattices have differences larger than one standard 
deviation, which is an indication of finite volume effects. Certainly, the inclusion of the 
errors due to the parameterization of the gluon dressing function used in the computation 
of Mpg^j(g^) improve the agreement between the different data sets. Note, however, that 
the figure [6] shows no clear systematic for the observed differences. Indeed, if the 64^ 
data is above the 80^ data, the 48^ data is below the largest lattice considered in the 
present work. 

Our interpretation of the results summarized in figure |6] being that the gluon can 
be viewed as having an effective mass for momenta below 1 GeV. For momenta above 1 
GeV, the gluon propagator follows closely the perturbative QCD prediction and becomes 
a massless boson. This behavior helps explaining why equation (jlj) could be used to 
described the lattice propagator for momenta as low as 1.5 GeV - see the renormalization 
procedure and results described in table [2] of section O 

3.2.2. Lattice M'^{q^) and Dyson-Schwinger Results: As discussed in |2], the 
nonperturbative solution of the QCD Dyson-Schwinger equations allows for a gluon 
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Table 5. Fits of the lattice data using equation ^T5\i (upper part of table) and ([TT]) 
(lower part of table) for mass definition. See text for details. 



dynamical generated mass. Such a gluon mass was investigated in detail within the so 
called decoupling solution of the DSE [H [101 HB [I2l US]. Typically, the DSE are solved 
numerically and the numerical solution is fitted to a given functional form for the full 
momentum range. By an appropriate choice of the fitting function, one can identify 
both the running gluon dressing function and mass. Here, we would like to check if the 
same functional forms are able to described the lattice propagator. Therefore, in this 
section, we will assume that the gluon propagator is given by the ([1]), with 

2\ ^0 



log — ^ 



7 ' 



(14) 



where 7 is the anomalous gluon dimension, mo a mass scale and r a numerical factor. 

For the effective gluon mass, we will consider two different functional forms. A 
more sophisticated expression, which incorporates the perturbative QCD ultraviolet 
prediction at high momenta [lOl [11], [121 [13] . 
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In the above expressions pi 
consider the expression 



1/2 and p2 = 5/2. Besides expression (ITS!) . we will also 



mn 



(17) 



for the effective gluon mass also used to described the numerical decoupling solution of 
the DSE. 

For the functional forms ( |T5l) and ( |T71) . the parameters are computed fitting the 
lattice data via minimization of the / d.o.f. For both expressions, we were not able to fit 
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Figure 7. Volume dependence of the parameters defined in ([T5|) and p7)) . 



all momenta range. This can be understood looking at the propagator data in ultraviolet 
region. Indeed, for large momenta, the statistical errors on the lattice propagator are 
quite small and the minimization process is constraint mainly by these points. Further, 
the renormalization process does not eliminate all the artifacts associated the finite 
lattice spacing. Indeed, a closer look at the propagator show that for g ~ 5 GeV and 
above, the propagators fluctuates slightly. Given the small statistical errors, this makes 
extremely difficult to fit the data unless the fluctuations are removed as, for example, 
in [251 [2ni |27]. Instead of smoothing our lattice data, we fitted ( !T5|) and (fTTll from 
g = up to Qmax, wherc q^ax was defined as the maximum momentum window where 
xVd.o.f. < 2. 

The fits are reported in table [51 Our first comment being that the fits are not able 
to distinguish between the two functional forms ( fT5i) and ( |T71) . The two functional forms 
describe quite well the lattice propagator in the nonperturbative region, i.e. from zero 
momentum up to ~ 4.2 GeV. Furthermore, the region which they are able to describe 
increases with the lattice volume. The exception being the smallest lattice, where the 
good fit for a relatively large Qmax is related to the smaller number of d.o.f., which makes 
the fit easier to perform. 

In the infrared region, (fTSll and (fTTll reduce to a massive type propagator whose 
mass is given by the mass parameter mo. Looking at the infrared mass computed 
in section 13.1.11 see table [3l it follows that the hard infrared effective gluon mass is 
slightly smaller than mg. However, mo is, within one standard deviation, compatible 
with the extrapolated infrared mass of 723(11) MeV. Looking at tuq obtained from the 
DSE equations, it turns out that their tuq is slightly smaller than the numbers reported 
here. For example, in [13] the authors get for mo = 612 MeV. Note, however, that our 
definition does not match exactly their definition. 
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Figure 8. Renornialized gluon propagator and q^D{q^) function from 80* lattice data 
and corresponding using for the mass squared the definition (jl7p . 

The volume dependence of the various parameters can be seen in figure [71 mo is well 
described by a linear function of 1/L, which can be used to extrapolate mo to infinite 
volume. The extrapolation to V — > oo gives mo = 671(9) MeV, for a /d.o.f. = 0.16. 
The extrapolated value is closer to the estimated infrared hard mass measured directly 
from the lattice data measured in section 13.1. II and to the DSE estimate for mo reported 
in [13]. For the remaining parameters, it is not clear which functional form should 
be used to extrapolate to infinite volume and, in the following, we quote the value 
computed from the largest lattice. Note that all the points for zq, A and r in figure [7] 
are compatible within two standard deviations. 

From the previous analysis one can conclude that the functional forms used to 
described the gluon propagator from the decoupling solution of the Dyson- Schwinger 
also describe the lattice data. Furthermore, they provide a definition for a dynamical 
generated gluon mass. 

The renormalized gluon propagator and corresponding q^D{q^) function for the 80^ 
lattice, together with the fits using ( ITTj) to parameterize the running gluon mass, are 
plotted in figure [H Figure [9] show the corresponding fitted dressing function and gluon 
mass. 

4. Results and Discussion 

In this work we have investigated if the Landau gauge lattice gluon propagator can 
be described by a massive type propagator. For the gluon mass itself, three different 
scenarios were considered: i) a constant infrared mass; ii) a constant ultraviolet mass; 
iii) a running mass in association with a running gluon dressing function. The mass was 
measured from the momentum space gluon propagator given by equation ([T]). 
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Figure 9. Gluon dressing function and gluon running mass given by (jl7p from the 
fits to the largest lattice volume. 

The interpretation of the infrared lattice gluon propagator as a massive propagator 
with constant M and Z was studied in section [3.1.11 Our results show that the lattice 
data is compatible with such a picture for momenta up to ~ 500 MeV, with an effective 
gluon mass around 650 - 700 MeV. If the zero momentum gluon propagator is included 
in the data sets, then the measured Z = 4.044(78) and M = 648(7) MeV. On the 
other hand, if one removes -D(O) from the analysis, Z = 4.49(10) and the effective gluon 
becomes slightly higher M = 723(11) MeV. 

The measured infrared hard mass reproduces the results reported in [5] for SU(3) 
simulations and is consistent with the quoted infrared mass value estimated for SU(2) 
in [28], where an M of 0.69(3) GeV or 0.68(4) GeV, depending if one includes or not 
include -D(O) on the analysis, was claimed. Furthermore, the gluon masses claimed 
above agree well with the estimate of the gluon mass from the gluon condensate (A^) 
obtained in [29], where the value 625(33) MeV was obtained, and is well within the 
interval of values estimated from phenomenology [1] . 

The interpretation of the lattice gluon propagator as massive type propagator with 
a constant mass for the ultraviolet region was checked in section I3.1.2[ It turns out that 
the lattice data cannot be fitted consistently by such a propagator. depends on 
the fitting range [qmin,qmax] and, for each Qmax, M'^ is not constant. This means that, 
in the ultraviolet, the lattice gluon propagator does not behave as a massive bosonic 
propagator with a non- vanishing constant mass. 

The case of a running gluon mass M^(g^) and a running dressing function Z{q^) 
was studied in section 13. 2[ Several definitions for M'^{q^) together with Z{q^) where 
investigated. In general, it turns out that is a decreasing function of q. Of course, 
the precise value for M'^{q^) depends on the chosen definition for Z{q^) and vice- versa. 
Here we considered two cases for the dressing function which allow for the computation 
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of the running mass: (i) a naive extension of the perturbative dressing function towards 
the infrared, keeping Z{q'^) always finite; (ii) the functional forms for M'^{q^) and Z^q^) 
used to fit the so-called decoupling type solution of the gluon-ghost Dyson- Schwinger 
equations. 

Our first computation of the gluon running mass considered a parameterization of 
the running dressing function which extended the perturbative Z towards the infrared 
region, while keeping Z{q^) finite for all momenta. The gluon mass defined in this way, 
Mpg^j(g^) is shown in figure El M(g^) runs from ~ 1.2 GeV for g = down to zero 
for q ~ 932 MeV. We would like to call the reader attention, that if instead of using 
the results reported in table H] for the larger fitting range, we used the results for the 
smallest fitting range, the values for Mperti.'f) would be reduced, starting below 1 GeV 
for the smallest momentum. However, the qualitative behavior will be similar to that 
observed in figure El 

After investigating Mp^^^^q"^), we have checked the compatibility between the 
functional forms used to described the decoupling type solution of the Dyson- Schwinger 
equations for the gluon propagator and the lattice propagators. These functional forms 
provide a parametrization of the gluon propagator and a definition for the running 
gluon mass. Moreover, they take into account the perturbative QCD predictions for 
the ultraviolet regime. The functional forms are able to describe quite well the lattice 
data over the entire non-perturbative regime. The lattice data can be fitted by the 
expressions considered in the present work from g = up to qmax, with q^ax increasing 
with the lattice volume. For our largest volume qmax = 4.13 GeV. For momenta above 
~ 4.13 GeV the remaining lattice spacing effects prevent a fit to the data. 

In what concerns the gluon mass, we found that it is well described by 



g2 + TTiQ ' 



with mo = 671(9) MeV. Our estimation of niQ agrees well with the estimation performed 
using the DSE and with is, within errors, compatible with the infrared hard mass when 
the -D(O) was included in the analysis. 

From the above results it follows that the interpretation of the gluon propagator 
as a massive type propagator with a momentum dependent mass and dressing function 
fits, quite well, the lattice QCD data. The gluon mass is a decreasing function of 
and becomes massless in the ultraviolet region. The nature of the gluon mass helps in 
the understanding of the remarkable mechanism of confinement, where the gluon mass 
may contribute to a Meissner-type effect in QCD. Moreover, given that one can define 
an effective gluon mass, it follows that the (A^) should be taken into account when 
investigating the non-perturbative djTiamics of QCD. 
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